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Abstract. Wg consider the Schrodinger equation with a random potential of 
the form 

Vix) = "'''^^^ - "''^^^ 
^ ' 4 2 

where ui is a Levy noise. We focus on the problem of computing the so-called 
complex Lyapunov exponent 

n := 7 - mN 

where A'' is the integrated density of states of the system, and 7 is the Lyapunov 
exponent. In the case where the Levy process is non-decreasing, we show that 
the calculation of O reduces to a Stieltjes moment problem, we ascertain the 
low-energy behaviour of the density of states in some generality, and relate 
it to the distributional properties of the Levy process. We review the known 
solvable cases — where Q can be expressed in terms of special functions — and 
discover a new one. 



1. Introduction 

The disordered model we consider is a particular case of the Schrodinger equation 

(1.1) + V{x)i: = Ell) . 

In this expression, if) is the unknown (wave) function of the independent variable 
a;, E is the energy parameter, and V is the random potential. If V has stationary 
increments, then the most accessible quantities describing the model's behaviour 
are the integrated density of states N{E), which counts the number of energy lev- 
els below E per unit length, and the Lyapunov exponent "f{E), whose reciprocal 
provides a measure of the localisation length at level E. These are conveniently 
expressed in terms of the limit 

lniP{x,E) 
iim 

where ip{-,E) is the particular solution of Equation (|l.ip satisfying ?/'(0, E) = and 
Tp'{Q,E) = 1. This limit is a self-averaging (non-random) quantity whose almost- 
sure value we denote by Q.{E). The relationship between the real numbers ^{E) 
and N{E) and the complex number Q,{E) is then (see [511 [551 [51] ) 

(1.2) Vl{E) ^-i{E) -iTTN{E) . 

We call f2 the complex Lyapunov exponent, or the characteristic function, of the 
disordered system (|l.ip . 
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Building on the pioneering work of Frisch & Lloyd [18], Kotani [26] made a 
rigorous and detailed study of the integrated density of states when the potential 
is a Levy noise. In the present paper, we shall instead examine the case where 

(1.3) Vi.) = 

and it is the superpotential w — rather than V itself — that is the Levy noise. Histor- 
ically, the study of the random supersymmetric case was initiated by Ovchinnikov 
& Erikhman [37] in the context of one-dimensional disordered semiconductors. The 
same model, expressed later on in terms of Equations (|1.15p and (|1.16p , may also 
be derived as the square of a Dirac operator with a random mass. As such, it is 
relevant in several contexts of condensed matter physics, including one-dimensional 
disordered metals, random spin chains and organic conductors; see [TSj [44] for re- 
cent reviews. From the point of view of Anderson localisation, the one-dimensional 
case is special, for localisation takes place as soon as there is any disorder [20l [27] , 
in contrast with the richer higher-dimensional case where a localisation transition 
can occur. Nevertheless, there are compelling reasons for undertaking this study: in 
the deterministic case, supersymmetry permits a detailed mathematical treatment 
of the spectral problem, leading in a few cases to exact results [15]. It is therefore 
natural to ask whether this favourable state of affairs persists in the presence of 
randomness. Recently, we showed how to extend the Frisch-Lloyd approach to 
a general class of random, point-like, scatterers, including one of supersymmetric 
type |14| , and one of our aims is to develop that preliminary work so as to bring out 
the very elegant structure peculiar to the supersymmetric case. Another powerful 
motivation is the equivalence between the supersymmetric model and the mathe- 
matical description of diffusion in a random environment [9j [T2j [29l [30l |4TJ [43] . In 
the remainder of this introduction, we review the terminology and give an outline 
of our main results. 

1.1. Levy processes. Any real- valued process with right-continuous, left- limited 
paths started at the origin, and stationary independent increments, is called a Levy 
process [1113]. Such a process, say W, is completely determined by its Levy exponent 
A, defined implicitly by 



(1.4) E (c'«^(-)) = e^^W . 

Furthermore, the Levy-Khintchine formula holds: 



(1.5) A{e) = 2.9 [iiJ0-6^] + 



- 1 - 



1 + r 



<dy) 



for some constants G M, 17 > 0, and some Levy measure to, i.e. a measure on M 
such that 

(1.6) / min{l, j/^} TO(dy) < cxj . 

Jr 

Conversely, given numbers iig and g, and a Levy measure to, there is a Levy process 
W with Levy exponent p.4p : fj.g, g and to are called the Levy characteristics of the 
process. 

Roughly speaking, the paths of a Levy process consist of intervals of drifted 
Brownian motion separated by jumps whose height and frequency are controlled by 
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the Levy measure. In the particular case where the Levy measure m is finite, i.e. 

(1.7) P ■~ "^(^2^) < ^ 

we call W an interlacing process; it may be expressed in the form 
(L8) W{x) ==2ax + 2y/gB{x) + J{x) 

where _B is a standard Brownian motion, the drift 2a is given by 

(1.9) 2a := 2/^5- / m{dy) 
and 

P{x) 

(1.10) Jix) = 

where P is a Poisson process of intensity p, and the hj are random variables with 
probability measure m/ p. The processes B and P, and the random variables hj, 
in these expressions are mutually independent. Figure [T] displays realisations of 
two important particular cases: (a) a standard Brownian motion, corresponding 
to the choice a = p = 0, and (b) a compound Poisson process with exponentially- 
distributed jumps hj, corresponding to the choice a = g = and p = 1. 




2 4 6 8 10 



(a) (b) 

Figure 1. W{x) against x for a particular realisation of (a) a 
standard Brownian motion and (b) a driftless compound Poisson 
process with exponentially-distributed jumps. 

Not every Levy measure is finite, and so interlacing processes do not exhaust 
the class of Levy processes; see Figure [2] for a particular realisation of a Levy 
process that is not interlacing, and [2] for a description of other important examples. 
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Nevertheless, every Levy process W may be approximated by an interlacing process 
We as follows: let e > 0, and define from m a new Levy measure by 

me(A) := m{Ae) where -.^ {y e A : \y\ > e} . 

The measure is finite. Hence the process with Levy characteristics ^g, g and 
Toe is an interlacing process, and converges to A pointwise as e — )■ 0. The sense 
in which the paths of approximate those of W can be made precise; see the proof 
of Theorem 1 in [4|. The intuitive picture suggested by this construction is that 
non-interlacing processes experience very small jumps with a very high frequency, 
as illustrated by Figure [5J 



6- 



4- 



2 




2 4 6 8 10 



Figure 2. W{x) against x for a particular realisation of the subor- 
dinator such that a = and the Levy measure is given by Equation 
([L26)) with (7 = 1. 

In what follows, particular attention will be paid to the subordinator case, where 
the Levy process is non-decreasing. For such processes, g = and the Levy measure 
does not charge (— oo, 0]. A subordinator is therefore of bounded variation on finite 
intervals and there holds 




min{l, y} TO(dy) < oo 



and 

/"OO 

(1.11) A{e)^2a+ (e'y" - 1) TO(d?/) 

Jo 

with a > 0. We shall see that this case, as in Kotani's study, affords a number of 
simplifications. 

1.2. The Frisch— Lloyd— Kotani study. Let us review the formal steps involved 
in the calculation of when the potential is the distributional derivative of a 
Levy process W. We refer the reader to the original papers [HI [26] for detailed 
explanations, hypotheses and proofs, and to [H] for heuristics. 



SUPERSYMMETRY AND DISORDER 



5 



Written in terms of the Riccati variable 

the Schrodingcr equation (jl.ip with the potential V = W = w becomes 

Z' ^-E- +w. 

Z is a Feller process, started at infinity, whose infinitesimal generator, say 
given by (see for instance [U [39] ) 



i^u) (z) = {2fig --E~ z^) iL'{z)+2giL"{z) + 



u{z + y) - u{z) - - 



yu'{z) 



V 



Z has a stationary distribution whose probability density, denoted /, is in the kernel 
of the adjoint '^^ . Hence 

{E~2iig + z^) !{z) + 2gf{z) + 



fit) At 



i + r 



m{dy) = N 



where the integrated density of states appears as a constant of integration; N may 
be determined by finding a positive integrable solution and imposing the normali- 
sation condition. The Lyapunov exponent 7 is then given by the Cauchy Principal 
Value integral 

/■oo 

1 = 4 zf{z)dz. 
<y ~Qo 

The density / contains more information than is strictly necessary to determine $7, 
and it is sometimes more convenient to work with its Fourier transform 



7(^)dz. 



/ satisfies the differential equation 

(1.12) f"{s)+c{is)f{s)^Ef{s), s^O, 

where 

A(is) 



(1.13) 



cis) - 



25[/i - s] 



1 



y 



It is then readily verified (see for instance 
(1.14) 



s 1 + 

Appendix B) that 



lidy) . 



il = 1— 



/(0+) 

where, with a slight abuse of notation, / now denotes any non-trivial solution of 
the differential equation (|1.12p such that 

lim f{s) = . 

This representation in terms of the decaying solution of a homogeneous second- 
order linear differential equation was used by Frisch & Lloyd as the basis of their 
numerical study of the density of states [18| . Explicit formulae were obtained by 
Halperin in the case of a Brownian motion |23| and by Nieuwenhuizen in the case of 
a compound Poisson process with a gamma distribution of the jumps [34]. Kotani 
performed a semiclassical analysis of the differential equation, in the limit of low 
energy, for the particular case where VF is a subordinator [55]. By using Langer's 
transformation, he was able to obtain an approximation of the decaying solution 
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valid uniformly in the variable s, from which the low-energy behaviour of N(E) 
could be deduced. 

1.3. The supersymmetric case. The Schrodinger equation with the supersym- 
metric potential ()1.3p can be recast as the Dirac system 

(1.15) -il;' - ^ ijj = Ve (f> 

(1.16) (j)' - ^ (I) = Ve ijj . 

The meaning of these equations when w — W' is a Levy noise becomes clear if we 
introduce an integrating factor: 



(1.17) 

ax 



W 

exp ( Y ) 



_Eexp 



2 



f W\ 1 r- ( ^\ 
exp ( — -\4> = V exp I — -\ 



Let us outline the modifications of the Frisch-Lloyd-Kotani approach required 
in the supersymmetric case. The appropriate Riccati variable is 

(1.19) Z:=-\/E-. 
Then 

(1.20) Z'^-E-Z^ + wZ 

and the Levy noise now appears as a multiplicative term; as we have written it, 
this stochastic equation should be interpreted in the sense of Stratonovich. In 
the supersymmetric case, the essential spectrum is contained in ]R+. It will be 
convenient to consider in the first instance the case E <Q,so that N{E) = and Z 
has a stationary distribution supported in whose density we denote again by /. 
To find an equation for /, we start from the simple observation that the logarithm 
of Z satisfies a stochastic equation in which the noise appears as an additive term, 
just as in the previous subsection. It is then straightforward to deduce 

(1.21) [E -2iigz + z^) f{z) + 2gz^ [z,f{z)] 
^ ' az 

/(t)dt + -^/(z)|m(dj/) = 7V. 

1.4. Exponentials of Levy processes. For E' = 0, Equation (|1.20p has a non- 
negative solution given by 

1 ^--W(x) 1 j r-\W(x)-W(t)\^^^ 




Z{x) Z(0) ^0 

In particular, if Z(0) = oo, we may use the stationarity of the increments of W to 
deduce that 



1 (law) 



e-'^Wdt. 



Z{x) JO 

So the reciprocal of the zero-energy Riccati variable has the law of an exponential 
of the Levy process W . The study of these exponentials has received a great 
deal of attention in the probability literature [5] , and its close connection with our 
disordered system makes it a rich source of ideas. In particular, Bertoin & Yor [6] 
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studied the moments of the exponentials and found for them a first-order recurrence 
relation which we proceed to extend to the case E 0. 

1.5. Main results and outline of the paper. In the supcrsymmctric case, be- 
cause the noise appears in Equation (jl.20p as a multiplicative term, the equation 
for the Fourier transform of / contains an awkward integral term. Following the 
example of Bertoin & Yor, it is in fact more expedient to work with the Mellin 
transform 

noc 

(1.22) f{s) / z'^f{z)dz. 

We derive in Sj2] the supersymmetric counterparts of the Frisch-Lloyd formulae 
(in2l) and (frT4| : for E < 0, 

(1-23) Ef{s + 1) - c{s)fis) + fis - 1) = 

and 

(1.24) n{E)^'M^EM. 

2 /(O) 

The coefficient c(s) appearing in the difference equation (|1.23|) was defined earlier 
by Equation ()1.13|) for s 7^ 0, and this definition may be extended to s = by taking 
the obvious limit. To find /(O) and /(I), one can either solve the integro-differential 
equation ()1.21|) for the distribution and then compute the relevant integrals or, 
alternatively, seek a positive solution of Equation (jl.23p satisfying the normalisation 
condition /(O) = 1. The problem of computing the density of states and the inverse 
localisation length in the "physical" region i? > is thus reduced to a problem of 
analytic continuation — in the energy variable — from ]R_ to the cut plane C\]R+. 
This formulation in terms of the Mellin transform sheds new light on the few solvable 
cases that have appeared in the literature; we review them in ^ 

The remaining sections are devoted to the subordinator case, where we have an 
effective tool for performing the analytic continuation. We show in 21 that 

(1.25) n{E) = 'M. 



2 - E 

c(l) 



c(2) + • • • 

for every complex E outside the essential spectrum, and we remark that the calcu- 
lation of the density of states is tantamount to solving a Stieltjes moment problem. 
In fj5l we look for instances where the general solution of the difference equation 
(|1.23|) is known explicitly. For the subordinator with Levy measure 

-qy 



(1.26) m{dy) = 



VI- 0-22/ _ 

we find, with the help of Masson |32| . that the complex Lyapunov exponent is 
expressible in terms of the parabolic cylinder functions. Then, in ^ we take up 
the problem of determining the low-energy behaviour of the integrated density of 
states and show how to compute the leading term in a semiclassical approximation; 
we also provide a list of the possible singular behaviours. Some of the implications 
of these results for diffusions in a random environment are discussed briefly. We 
end the paper in Sj7]with a few concluding remarks. 
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2. Characteristic function and Mellin transform 



This section provides a derivation of Equations (jl.23|) and (|1.24p . For the partic- 
ular ease where is a subordiuator with finite means and a finite Levy measure, 
this derivation may be made completely rigorous by adapting the arguments of 
Frisch & Lloyd [l^ and Kotani [2^. In the general case, however, we shall be 
content to view the formulae as merely plausible. 

2.1. The diflference equation. First, as a consequence of Sato's Theorem 25.3 
[40j . we have the following criterion for the existence of the coefficient c(s) defined 
by Equation (|1.13p : for s 7^ 0, c(s) exists if and only if 



e m{dy) < 00 . 



\y\>i 

Also, it is readily seen by direct calculation that c(0) exists if and only if 
(2.1) J ym{dy) < 00 

and, in case of existence, 

c(0) =E(VK(1)) . 

In particular, c(s) exists for every s > if is a subordiuator with finite means. 

Next, to derive the difference equation, multiply Equation (jl.2ip by and 
integrate over z. We use 

z— [zfiz)] dz = z'-'zf{z) -(1 ~ s) / Z-' zj{z) dz = (s - l)/(s - 1) 

where we have assumed that 



(2.2) 
Also, 



lim z^^'fiz) = and lim z^^"/(z) = . 

z— >0+ 2— >oo 
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fit) dt 



yz 



dz 



f{t)dt 



V 



1 - s 



l + y2 
f{t)dt 



1 



1 - ,s 
1 



z'-^ [f{ze-y) - f{z)] dz + -^f{s - 1) 



1 



z^-'f{zc-y)c-ydz- 



1 



1 



/(«-!) + 



By making the obvious substitution in the integral, we obtain 

yz 















Jo 


J z 



fit) dt 
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When we put these resuhs together, we eventually find that /(s) solves the differ- 
ence equation (|1.23p . 

Let us now look back on Condition (|2.2I) in the case s > 0. The fact that the 
limit at infinity vanishes is a consequence of the Rice formula 

lim z'^f{z) = N = 0. 

z— f oo 

That the limit at zero vanishes is obvious if ly is a subordinator and E = —k"^ < 0, 
because the support of f{z) is contained in [k, oo). 



2.2. The Lyapunov exponent. For E ^ Q, 



— In — In 

X X 



Z{x) 



1 



1 



1 



In J\E\ + - In - - In \Z{x)\ . 

XXX 

If the Levy process is not a pure drift then the Riccati process will behave ergodi- 
cally, and so the alternative formula 

1 



(2.3) 



liE) 



lim — In 10(0;) I 

i-oo X 



will hold for E ^ 0. Our aim is to use this formula to express the Lyapunov 
exponent in terms of the Mellin transform of the Riccati variable. 

We shall suppose that W is an interlacing process of the form (|1.8m.l0]) , with a 
finite mean at a; = 1. By Equation (j2.1[) this implies that the jumps hj themselves 
have a finite mean. The Dirac system 15111. 16|) is then 



^ 00 



Hi}}. 



Here 



xq := and Xj+i := Xj + (.j 
where the Ij are independent and exponentially distributed with parameter p. For 

X G {xj, Xj+i), we have 

Hence 



ln|0(a;j+i-)| =ln|0(a;,+)| + 



^3+1 



dy 



m 

In \(l){xj+)\ + alj + ^ [B{xj+i] 



B(x,)] 



^E 



z{y) 



dy. 



On the other hand, by making use of Equation (jl.lSp . we find 

' h-j 



exp 



and it follows that 



In = In +a^, + [B{x,+i) - B{xj)] + ^ + 



Xj + l 



-E 



dy. 
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Then, by summing over j, we obtain 

fP{x) . 
ax + ^B{x) + lf^hA+^J^ 



In the hmit as x ^ oo, this becomes 

(2.4) + r^/(.)d.. 

This estabUshes the formula (|1.24p for the complex Lyapunov exponent when W 
is an interlacing process with finite means. We expect this formula to hold more 
generally even in cases where the Levy measure is not finite since, as mentioned in 
m.ll every such measure may be approximated by a finite measure. 

3. Some known solvable cases 

There are very few known cases where the complex Lyapunov exponent may be 
expressed in terms of familiar functions. Whereas the examples we review here 
were discovered by solving the integro-differential equation (|1.2ip directly, our pre- 
sentation will emphasise the alternative approach based on solving the difference 
equation satisfied by the Mellin transform. 

Example 1. When the Levy process is a pure drift, the process is deterministic 
and we have 

A(e') = 2106* and c(s) = 2a. 
The difference equation (|1.23p takes the simple form 

(3.1) Ef{s + 1) - 2af{s) + f{s - 1) = . 

Its general solution is 



c+ [a + - e) + c_ (a - - E 



The Mellin transform of the invariant density is positive for every s and equals 
unity for s = 0. Hence 

and, by using Equation (|2.4p . we find 



(3.2) n = - E . 

The spectrum is the interval E > a^; there, we have 

N{E) = -VE-a^ and j{E) = . 

TT 

Hence 7 vanishes in the spectrum — reflecting the fact that, for the non random 
potential V{x) = a^, there is no localisation. 

Example 2. The case where the Levy process is a Brownian motion with drift has 
been studied independently by Ovchinnikov and Erikhman |37| and Bouchaud et 
al. [9]. In this case, 

A{9) = 2g [ifiO - e^] and c(s) 2g - s] . 

Hence /(s) satisfies the difference equation 

Efis + l) + 2g[s- ^i] f{s) + fis - 1) = . 
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We recognise the reccurence relation satisfied by the modified Bessel functions. The 
general solution is 



and we obtain a positive solution by taking c/ = 0. 

Alternatively, we can solve the equation (|1.2ip for the invariant density. Since, in 
this case, the Levy measure is identically zero, this equation reduces to a differential 
equation. The relevant solution is (before normalisation) 



f{z) = ^ cxp 



1 

2^ 



E 



and its Mellin transform 



^-E 



L(0,< 



^/g 



agrees, as it should, with the expression found earlier by solving the difference 
equation. Hence, by Formula (|1.22|) . 



ri(-fc^) = ^ig + k 



Ki-^ jk/g) 
K-f. (k/g) 



where we have set E = — fc^. The analytic continuation of O from 
be done "by hand" and consists of replacing k by — ifc: 



to may 



n(k^ + iO- 



fj-g 



H\'},{k/g) 
^ H^'lik/g) 



Expressions for TV and 7 follow easily. In particular, for /i = (the driftless case), 
we find 



NiE) 



2g 

\nE 



and ~ — — In £■ as 

l{E) 2g 



0- 



So the density of states has a so-called Dyson-type singularity at the bottom edge 
of the spectrum, and the model exhibits a very interesting "delocalisation" phe- 
nomenon there [HI [131 SI] • 



Example 3. Let the Levy exponent be 



/ [c'"'' - 1] qc-'^y dy 
Jo 



i9 



P 



q-ie 



This corresponds to a driftless compound Poisson process of intensity p where the 
jumps are exponentially distributed with parameter q. Then 



c{s) = 



q + s 



s > 0. 



Set 

Masson [ 
given by 



E=~k^, k>0. 
showed that the general solution of the difference equation p.23|) is 

C+W+(s) + C-U-{s) 
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where 

u±(s) :=(±fc)' 



X2Fi(q + s + l,± 



P 



1± 



P 



-1 



'2k ' ' " ■ - 2fc' 
and B is the Euler beta function. To obtain a positive solution, we set c_ = 0. The 
formula cqrefmcUinForniula for the complex Lyapunov exponent then agrees with 
the calculation of our previous article [13] , where we have worked directly from the 
invariant density (before normalisation), namely 



-k 



- fc2 \z+k 



)(^)- 



After normalisation we obtain the Lyapunov exponent from Equation (|2.4p [14j . 
For negative energies 7 coincides with the characteristic function H. After analytic 
continuation we find the integrated density of states. The low-energy behaviour is 
(see Bienaime [7]) 

(3.3) N{E) - ^"^ , E-'i exp 



and 



7 ' 



P_ 

2q 



as E^O- 



Remark 3.1. The low-energy behaviour of the density of states in the supersym- 
metric case is quite different from the usual Lifshits behaviour I28| . namely 



N{E) - exp 



as ^ 0- 



The physical reasons for this difference will be explained in a forthcoming paper. 



4. Continued fractions 

There is a well-known connection between linear second-order difference equa- 
tions and continued fractions [THl ; we exploit it to derive a continued fraction 
that may, at least in the subordinator case, be used to compute the complex Lya- 
punov exponent. 

4.1. Pincherle's Theorem. We say that a non-zero solution Un of the difference 
equation 

(4.1) w„+i = a„u„_i + &„u„ , n e N , 

is recessive (or minimal or subdominant) if, for every other linearly independent 
solution w„, there holds 

Un = o (vn) as n ~> CXJ . 
Pincherle's Theorem says that the continued fraction 



Ol 



62 + • • • 

converges if and only if the difference equation (|4.ip has a recessive solution u„ with 
uq 7^ 0. In case of convergence, the limit is —ui/uq where m„ is any such recessive 
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solution. The proof of this weh-known resuh is elementary, but we include it here 
because it contains information that will be useful later on. Introduce two particular 
solutions pn and g„ of the difference equation (|4.ip satisfying 

Po = qi = 1 and qo = pi = . 

These solutions are obviously linearly independent, and every other solution may 
be expressed as a linear combination of them. Set 

Kn:=—, n€N. 

q-a 

To say that the continued fraction converges is to say that has a (finite) limit 
as n oo. First, suppose that the continued fraction converges; call its limit K^a 
and set 

■ — qn -^oo Pn • 

It then follows from this definition of the sequence m„ that mq = — 1 and 

Aoo = • 

Uq 

Furthermore, m„ solves the difference equation. To show that u„ is recessive, sup- 
pose that Vn is any other linearly independent solution. We can express w„ in the 
form 

Vn = au„ + Pqn 
for some number a and some non-zero number /3. Then 

— = q: + /3 — = a + — ^ > /3 X oo . 

Un Un Aoo — ^ n^oo 

This shows that u„ is recessive. Conversely, suppose that the difference equation 
has a recessive solution u„ such that uq ^ 0. Since go = 0, the particular solutions 
qn and u„ are linearly independent. Therefore 

hm ^ = . 

ri-i-oo qn 

On the other hand, we can express u„ as a linear combination of p„ and g„: 
It follows that 

Vn 1 / Un \ Ui 
— = — Ul > . 

qn Uq \qn J "^oo Uq 

To complete the proof, there only remains to observe that there cannot be two 
linearly independent solutions that are recessive. 

Remark 4.1. The proof shows in particular that, if K converges to K^o, then every 
recessive solution of the difference equation is a multiple of 



qn ^oo Pn 
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4.2. The subordinator case. Let us now elaborate the relevance of Pincherle's 
Theorem to our study. For convenience, we shall omit the drift, so that the Levy 
exponent is given by Equation (jl.lip with a = 0, and the coefficient c(n) by 

f°° 1 — e^"^ 

(4.2) c(n) = / m(dy) . 

Jo n 

Set 

(4.3) Un-.^ E''f{n), neN. 

The fact that /(s) solves the difference equation (|1.23p for s € N implies that Un 
solves 

(4.4) Un+l = ~EUn-l + c{n)un , nSN. 
The corresponding continued fraction is 



(4.5) K{E) := 



-E 



c(l) 



-E 



- E 

c(2) 



c(3) + • • • 

We now observe that the coefficients c{n) are obviously positive. This means 
that the continued fraction (|4.5p is of the kind studied by Stieltjes [31 IMl HH US . 
Stieltjes showed that a necessary and sufficient condition for the convergence of the 
continued fraction in C\R+ is that the series 

oo 

^c(n) 



diverge. Now, by Equation (j4.2p . 

1 — c^"^ f°° 1 — c^"y 
c("-) — / m(dy) + / m{dy) 



> I ym{dy)+ / m{dy) 



n 

oo 

min{l, y} m{dy) as n — > oo . 



Hence the continued fraction converges for every E G C\]R+, and the limit Kao is 
a function of E analytic in C\R+ . 

Next, let us show that E'^ f{n) is a recessive solution of the difference equation 
when E <Q: since it is a solution, we may write 



E''f{n) = a [qnKooiE] - Pn] + Pqu 

for some constants a and /3. Note in particular that the g„ are positive, because 
the c{n) are positive. Suppose that /3 7^ 0. Then, since qnK^ — Pn is recessive, we 
must have 

E^ f{n) ^ Pqn as n — > 00 . 
But this is absurd because the right-hand side is of the same sign for every n 
whereas, since E <Q and /(n) is by construction positive, the left-hand side alter- 
nates in sign. Hence /3 = and the recessiveness of E"f{ri) follows from Remark [4.1l 
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We deduce from Pincherle's Theorem that 



for E' < 0, and thence for every E e 



by analytic continuation. 



Remark 4.2. This representation of the complex Lyapunov exponent in terms of the 
continued fraction associated with the difference equation for the Mellin transform 
is not valid for every Levy process. For the Brownian motion with drift of Example 
m f{s) is not recessive, and the continued fraction converges to a ratio of I^. 

Remark 4.3. Continued fractions have been used before in the study of other 
one-dimensional disordered systems — most notably by Nieuwenhuizen [SI] and 
Nieuwenhuizen & Luck |35j . In these works, however, the continued fraction arose 
only in very special cases where the noise or randomness could be modelled in terms 
of the exponential and other closely related distributions. 

4.3. The density of states in the subordinator case. In particular, we have 
the formula 



(4.6) 



N{E) 



1 



Im 



Ml 

uo 



1 



lm[K^{E + iO+)] , E>0. 



On the other hand, since the c(n) are all positive, the limit of the continued fraction 
is the Stieltjes transform of a certain measure, say k, on the half-line M+ and we 
may write (see for instance [31 ) 



(4.7) 



E 



K{dt) 

t-E 



E e 



There is a simple relationship between the Stieltjes measure n and the integrated 
density of states TV: by the Stieltjes-Perron inversion formula, for every open in- 
terval (a, b) C R+ 



1 r 

K{a, &) = / Im 

We then deduce from Equation (|4.6p 
(4.8) 



dt. 



, . , dt 



Now, as is well-known, the coefficients c(n) appearing in the continued fraction K 
may be expressed in terms of the positive integer moments of the measure k; see |42| 
for detailed formulae. So the problem of finding the density of states is essentially 
equivalent to a Stieltjes moment problem. 



5. A NEW SOLVABLE CASE 



Having reviewed in ^ the few examples where the complex Lyapunov exponent 
is known explicitly, we now look for new cases. 
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5.1. Masson's difference equation. Masson's catalog of the solutions of the dif- 
ference equation 



(5.1) 



Vn+i - (z - Dn)vn - {Am?- + Bn + C)ij„_i = , n e N , 



expressed in terms of the hypcrgcomctric function or its confluent limits, is particu- 
larly valuable in this respect [32] . The correspondence between Masson's difference 
equation and (|4.4[) is as follows: set 



c(n) 



and 
(5.2) 



Then u„ solves the difference equation (|4.4p if and only if Vn solves the difference 
equation 

1 1 

Vn+l P=^^n '"ri-l = . 



C„_iC„ 

This is Masson's difference equation ()5.ip with D — Q and z ~ 1/^ —E provided 
that 



(5.3) 



1 



-ic„ 



An? + Bn + C 



There are two possibilities: either 

c{n) r 



Ai 



or c{n) 



'En 



The first case leads to solutions expressible in terms of the hypcrgcomctric function; 
this is our earlier Example [S] The second case, developed in the next subsection, 
leads to solutions expressible in terms of parabolic cylinder functions; we follow 
Masson's terminology and call it the "Hermite case" . 



5.2. The Hermite case. Let the Levy process be the subordinator with Levy 
measure 



(5.4) 



m(dy) 



2p d 
Vtt dy 



Vi - c-2y 



dy , g > . 



The only dimensional parameter in this expression is p, which has the dimension of 
an inverse length, or equivalently of \fE. For convenience, we shall for the present 
set it to unity, and reintroduce it later whenever it becomes relevant. To work out 
the corresponding Levy exponent, introduce the tail of the Levy measure: 



m{dt) 



,y Vl - e-2y 

This makes it clear that m is not a finite measure and that 

dy 



(5.5) 



m{dy) 



y 



3/2 



as y ^ 
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We have 

m = 

Jo 
2i0 



(c'"^ - 1) m{dy) ^lO Ay c'"^ ?7i(y, oo) 



Jo 



' dt . = — =B ' 



VI - e-2i/ V^7o v^Frt2 V^ V 2 '2 

Here, we have made use of Formula 1 in §8.380 of [51]. Hence 

r 



A(6') = \e- 



q—lt 
~2~ 



r 



q-ifl+l 
2 



We note for future reference the large behaviour 
(5.6) A{e) V-2i6l as 61 ^ oo 

which is directly related to the small y behaviour of the Levy measure in Equation 
|). The coefficients appearing in the difference equation (|4.4p are given by 



(5.7) 
Hence 



c(n) = c„ := 



n + q - 1 



Cn—lC-n 2 

Let u„ and Vn be related by Equation (|5.2I) . Then Equation (|5.3p says that 
u„ solves the difference equation (j4.4p if and only if u„ solves Masson's difference 
equation with 

z = -^=, A = 0, B^l, C = ^— i and Z) = . 

2' 2 

Masson's Theorem 4 [31] then says that the difference equation (|4.4p has two linearly 
independent solutions given by 

(5.8) - ' + 



This fact is easily verified by using the recurrence relation satisfied by the parabolic 
cylinder functions (see [T7], Vol. 2, Chapter VHI): 

Dy+i (x) - xD,, {x) + vD^^i{x) = . 

Furthermore, is recessive. By using the identity 



1 



V{2z)^-^2'^-'T{z)T\^z + - 
we deduce the following expression for the Mellin transform of the invariant density: 



1 



(5.9) 



r(^±f±i ) D- 



5-9 



^2/i; 



Therefore, after re-introducing the parameter p, we obtain 



(5 10) o(i,)-^ ^(i) I ^(i) 

^ ^ ^ "^'^^-2r(^)+ r(^) 
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Using D-g{0) = 2-«/20F/r(i±2) wc obtain 

(5.11) !](£;)- 7oo + as£;^-oo. 

The constant term is given by 



(5.12) 
where 



1 2^ r(^)^-§r(§)^ 



:(0) =E(VK(1)) =p 



An analytic continuation to positive energies shows that 

The Stieltjes measure of a ratio of parabohc cylinder functions is given exphcitly 
in Ismail & Kelker [25], Theorem 1.4. That result, combined with Equation (|4.8p . 
leads to the following explicit formula for the integrated density of states: 



(5.13) 



N{E) 



r(§) 



20Fr(g)r(^ 



The low energy behaviour is 



(5.14) 



N{E) 



r,29 



q yip^2/E 



Plots of N and 7 for various values of the parameter q are shown in Figure |3l 

Remark 5.1. Let E' < and set k ~ \J —E. By using Equation (7.727) in [21], we 
can invert the Mellin transform / to obtain 



(5.15) 
where 



/(z)-C(fc,<z)- 



(z2-/c2)^ 



exp 



l(fc,oo)(2) 



C(fc,q) :=2 



Remark 5.2. The limit g — > 0+ is singular, since 



E(VF(1)) ==c(0) ==p 



r(f) 



00 as g — >■ 0+ . 



On the other hand, if we set g = in Equation (|5.4p . we obtain the Levy measure 
of the process with exponent 



This process is well-behaved; indeed 

c{n) = p 



r(t) 



r (ii±i) ^r, 
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and so 

E(M^(1)) = c(0) 4=21n2 for g = . 

Unfortunately, Equation (|5.3p docs not hold in this case, and so wc cannot use 
Masson's results. 




Figure 3. Integrated density of states and Lyapunov exponent for 
the Hermite case with q ~ 0.2 (black continuous line), 0.5 (green 
dotted line), 1 (blue dashed), 2 (red dashed- dotted) and 5 (black 
continuous) . 



Remark 5.3. One reason for the importance of the subordinator subclass of the 
Levy processes is that they arise naturally in the study of diffusion processes. More 
precisely, let X = X(t)he a. diffusion process on R, started at the origin, and driven 
by the Langevin equation 

X = v{X) + r] 

where 77 is white noise and v is the deterministic drift function. Then the inverse 
local time of X is a process, say Tt, that measures how much time elapses before X 
has spent a total time t in the vicinity of its starting point. It is a well-known fact 
that r is a subordinator. Given the drift v, one may in principle compute the Levy 
exponent of the inverse local time by studying a certain Sturm-Liouville problem 
associated with the diffusion. In particular, for the drift 

v{x) ~ —X , 

which corresponds to the familiar Ornstein-Uhlenbeck diffusion process, one finds 
that the Levy measure of t is exactly the measure m in Equation (|5.4p for q = 
1; see for instance It is unclear whether this intriguing coincidence is the 

manifestation of a deeper connection between the Ornstein-Uhlenbeck process and 
the disordered model. 

6. The low-energy behaviour of the density of states 

Throughout this section, we suppose that the Levy process W is a subordina- 
tor. Whereas Kotani's study rested on a semiclassical analysis of a differential 
equation, we must, in the supersymmetric case, work with a difference equation. 
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Finite-difference versions of the WKB metliod were expounded by Dingle & Mor- 
gan |16| and Braun and we sliall draw on tlieir work in what follows. It will 
be convenient to recast the difference equation (j4.4l) in the form 



(6.1) 

obtained by setting 

The basic formula 
(6.2) 



Un+l + Un-l 



c(n) 



Un = E^y^. 
for the integrated density of states becomes 



N{E) 



E . 



Im 



yi{E + iG+) 



where y„ {E) is any recessive solution of Equation (|6.ip . 



6.1. The remainders. The starting point of our analysis is the represention of 
the recessive solution in terms of the remainders 

E 

(6.3) z„ := 



-E 

c{n) + — ^ 

c(n + 1) + rr 

^ ' c(n + 2) + • • • 

These remainders are certainly well-defined numbers for E ^ M+ because the con- 
tinued fraction on the right-hand side converges. We note the obvious identity 

E 

(6.4) Zr, = -J- . 

Set 

n 

If we take the empty product to equal 1, it is easy to show by induction on n that 

E^y^ ~ E^^f{n) for every n>0. 

This shows that y+ is a recessive solution of the difference equation (|6.ip . We 
remark also that, since c{n) decays to zero, z„ and Zn+i have a common limit, 
namely —\/—E, as n — >■ cx). Hence, replacing Zn+i by z„ in Equation ()6.4|) yields 



cjn) - ^cjnf - iE 
(o.d) z„ ^ as n — oo . 

We select this particular solution of the quadratic equation because it yields the 
correct behaviour of z„ in the limit i? — )- 0; see Equation (|6.4p . 

Henceforth, we suppose that E > and set k = \fE. The expression y^iE) 
should then be interpreted as y^iE + iO+). 
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6.2. An alternative formula for the density of states. The difference equation 
()6.ip has a turning point Uk, defined to be the largest integer not exceeding the root 
of 

c(n) = 2k. 

We proceed to express the integrated density of states in terms of the recessive 
solution at the turning point. 

Let < E < c(l). For n > Uk, we have 



where 



cifi) 

On = arccos — — — > . 
2k 



It then follows from Equation (|6.6p that, for n > Uk and E > 0, there holds 



(6.7) 



yt - Vn, exp 

j>nk 



as n — > oo . 



Since c(n) is real, the complex conjugate 



Vn ■= yt 



is also a solution of the difference equation (|6.ip . Furthermore, and are 
linearly independent. Introduce their Wronskian 



Vn-l Vn 

yt 



Then 



N{E) = — Im 
27r 



2ti 



ImWi 



On the one hand, by making use of the difference equation (|6.1|) . it is easy to see 
that the Wronskian of any two linearly independent solutions satisfies 

Wn = Wl 

for every n G N. On the other hand, by construction, has the asymptotic 
behaviour (|6.7p . Hence 

-\ynS [e''" - = I y+ I ' [2i sin 0„] ^ | y+ | ' 2i as n ^ oo . 

We deduce 

Wi = |'2i 

and so 

(6.8) iV(£;) = ^|y+r. 
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6.3. Semiclassical approximation. There remains to find the small k behaviour 
of y^^. From this point on, we shall not aim at a completely rigorous treatment, 
but rather at a simple informative estimate based on heuristic considerations. 
Following Braun we introduce a small parameter h > such that 

h ^ a and only if fc — > . 

The relationship between h and E is assumed to be such that there exist functions 
C and Z defined on M+ satisfying 

(6.9) C{nh) = ^ and Z{nh) = ^ . 

Detailed examples will be given in due course. Equation ()6.4p for the remainder 
may then be expressed as 

1 



(6.10) 2C{x)-Zix + h) = 

Z{x) 

where 

X = nh . 

We look for a solution of Equation (|6.10p of the form 

(6.11) Z{x) =exp[Soix) +0{h)] as h ^ . 

When we substitute into Equation (|6.10p , use the MacLaurin series of the exponen- 
tial function, and equate like powers of h, we find, for the leading term, 

p-^O _|_ a — So 

(6.12) Y 

Then, by making use of the Euler-MacLaurin formula, we obtain 

(6.13) ?/„+ = n^(j-M=exp{i |%o(t)di+^5^M±M^ + o(i)| . 

In order to have 

Z{nh) ^ ——— as fc — > 
c(n) 

for n < Uk, we must select the following solution of Equation (|6.12p : 

So{x) = — arccosh [C(a:)] . 

Then 

= exp 1"-^ y arccosh [C{t)] dt - -arccosh [C{h)] + 0(1) 
where Xk is defined by 

(6.14) Cixk) = 1 . 

In view of Formula (|6.8|) . this leads to the estimate 

(6.15) In = -- / arccosh [C(t)] dt - arccosh [C(/i)] + Infc + 0(1) 

" Jh 

as k — > 0. The validity of this estimate is open to question, for we have assumed 
implicitly that the 0(1) correction term in Equation ()6.13|) is bounded uniformly 
for h < X < Xk- We proceed to discuss some examples where this assumption seems 
justified. 
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Example 4. For the Levy process of Example |3l set h = k. Then 

c{n) p 



2k 2 {nh + qh) 



C{nh) . 



The undesirable dependence of the function C on h could be removed by shifting the 
index n. As it is, however, Formula (j6.15p makes sense and may be used directly; 
shifting the index would eventually lead to the same result. We have 



/ arccosh 


P 


2{t + qh)_ 


Ih 



y=c{t) 

4- 



arccosh y 

y 



arctan ■ 



We deduce 



-^|-2(g+l) Inh + 0{1) 



a.s E ^0+ 



N{E) - AE-'' exp 

for some constant A independent of E. This agrees with Equation p.3p . 
Remark 6.1. We remark that, for every finite Levy measure, 



1 „ Q-^iy 1 

m{dy) — j m{dy) as n — cx3 . 



c(n) := 

Therefore the asymptotic relation 

TT 1 f°° 

(6.16) \nN{E) ^ m{dy) a.s E 

2 y/E Jo 



0- 



holds more generally for every subordinator with a finite Levy measure. 
Example 5. For the solvable model of we have 



c{n) 



Strictly speaking, the corresponding function C is defined implicitly by Equation 
(|6.9p . but a more tractable expression may be obtained by using the identity 

o 

c(n — 1) c(n) 



Then 



C{nhf 



1 



2fc2 In + q 



n + q — 1 
- + 0{h^) as /i 



So, for the purpose of our calculation, we may take 
h 2fc2 and C{x) := 



1 



h{q-^) 
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Then 



arccosh [C{t)] dt = 



V=C{t) ( 

-I 4 I ^yy^ - 1 1 

■' — — ^arccosny 



2y 2y2 



2 

h 



9 + - ) ln/i + 0(l) as /i^ 0. 



Formula (|6.15|) then yields 

N{E) ^ AEi~''e~i as E ^ 0+ . 
This is in agreement with the exact result found in see Equation (j5.14p . 

Example 6. The so-called alpha-stable suhordinator has for its Levy measure 

a dy 



m{dy) = p ■ 



< a < 1 , 



r(l - a) yi+" 
The corresponding Levy exponent is 

(6.17) A{e) = -p{-ie)" . 

Some of the properties of alpha-stable subordinators are reviewed in the appendix. 
For p = 1, the continued fraction coefficients are 

c{n) = 

and so we set 

h := (2fc)T^ and C{x) := x""^ . 

Then 

2 r^'' 2 

arccosh"^ [C(t)] dt = --^a - 2(1 - a) In /i + 0(1) as ft. ^ 



where we have set 

(1 - a)2-Fi( ^, 



r;i + 



1 



^2'2(l-a)' 2(l-a)' 
We deduce from Formula (|6.15p 



1 = 



2 ' 2(1-q) ' 2 



3.18) N{E) ^ AE^ exp 



2T^ 



as E^O^ 



6.4. The case of infinite Levy measure. Let us now discuss the case 




in greater generality. A first interesting observation is that the two instances en- 
countered earlier, namely the Hermite case and the alpha-stable case with a = 1/2, 
have led to the same leading behaviour of the integrated density of states at low en- 
ergy; compare Equations (|5.14p and (|6.18|) . The mathematical explanation is that 
the corresponding Levy measures have the same singularity at y = 0; the fact that 
their tails are very different has little bearing on the low-energy behaviour of N. 
The physical interpretation is that, for an infinite Levy measure m with support 
R+, the low energy eigenstates of the supersymmetric Hamiltonian are strongly 
affected by the frequent small jumps but depend weakly on the rare large jumps 
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of the Levy process. This contrasts with the case of a finite measure, where the 
exponent appearing in Formula (|6.16p indicates that there is no such distinction. 

Turning then to the analysis of the low-cncrgy states of the supcrsymmetric 
Hamiltonian, consider the finite interval [0,£]. In the absence of any potential, this 
interval supports an eigenstatc of (kinetic) energy i^kin ~ l/i'^- A low-energy state 
E w l/£^ is possible if, in the interval < x < £, the potential V = — w'/2 

remains sufficiently small to ensure that V <^ -Ekin- Because W{x) is a subordinator, 
we expect that the potential's average value is roughly proportional to the square of 
W{£)/£. Hence, if W{£) does not exceed some small value, say e, then the potential 
brings in a contribution V sa {£/£)'^ to the energy. Equating £ and 1/\/E, wc deduce 

(6.19) \nN{E) ^InF (w (i/Ve) < e] as E ^ 0+ . 



We argued earlier that, in this limit, only the frequent small jumps of the process 
matter. With this in mind, let us begin with the case of a Levy noise whose measure 
exhibits the "alpha-stable singularity" 

dy 

m{dy) p , < a < 1 , as y ^ 0+. 

In the appendix, we describe the implications of this singularity for the form of the 
probability density function, say fw{' ! 2;), of the random variable W{x); we find 

Equation (|A.4p in the appendix, which describes the small u behaviour of this 
density, then leads to 

(6.20) liiN{E) r^ln fw (u;1/Ve^ du E~^i^ as E ^ 0+ . 

Modulo a constant factor on the right-hand side, these heuristic arguments extend 
our earlier result (|6.18p . 

Finally, let us mention very briefly the limiting case a — )■ 0, where the Levy 
measure satisfies 

(6.21) m(dy) - — as y ^ 0+ . 

y 

The corresponding process W shares this singularity with the so-called gamma 
subordinator; see the appendix. By the same heuristic arguments as before, we 
arrive once again at Equation (|6.19p . and the result is 

(6.22) lniV(i;) ^ -Mi/^ asy^O+. 

Ve 



6.5. Comparison with the Kotani case V = w. It is interesting to compare 
these results with those obtained by Kotani [26] for the random Schrodinger equa- 
tion 

+ wiIj = Ei/j. 

For the compound Poisson process of Example HI Kotani showed rigorously that 



3.23) N{E)r^Aeyi-p 



7e 



as E ^ 0, 



26 



ALAIN COMTET, CHRISTOPHE TEXIER, AND YVES TOURIGNY 



where p is defined as before by (|1.7p . Tliis is tlie weil-fcnown Lifshits singularity 
[28] . For the alpha-stable subordinator of Example El he found 



(6.24) 
where 



N{E) ~ Aexp 



asE^O 



r(l-a) 



B 



\^ 1-a ' 2 



Kotani also obtained exact expressions for the constant factor A; see [55] , Theorem 
4.7. In the supersymmetric case, the constant factor may in principle be obtained 
by including the next term in the WKB expansion of the remainder Z, but the 
calculations involved are tedious. 

The asymptotic behaviour (|6.24p may be recovered by a heuristic argument sim- 
ilar to the one developed in the supersymmetric case. The main difference is that, 
whereas in the supersymmetric case V = w'^ /A—w' /2 the Levy process is dimension- 
less, in the case V = w it has the dimension of the reciprocal of length. Consider 
as before an interval [0,i], and suppose that the Levy process does not increase 
beyond a certain threshold e in that interval. The contribution of the potential 



energy now reads V ~ e/i. Therefore a state of low energy E w iJkin 
possible if y ^ Skin i-C. e <C 1/^ w a/B- Hence 



is 



In N{E) - In I 



(w{i/^ 



Ve] < e 



as E^O- 



The different exponent, compared to Equation (|6.20p . comes from the fact that the 
threshold e is now dimensional and scales with energy as \/E. The case m{dy) ^ 
dy/y for y — > 0+ leads to the same behaviour as in the supersymmetric case; see 
Equation ([QS]) . 

6.6. Classical diffusion. The formal equivalence between supersymmetric quan- 
tum mechanics and classical diffusion is well-known and has proved fruitful in the 
study of diffusion in a random environment [HI [301 US] ■ The purpose of this 
subsection is to outline the most immediate implications of our results for this field 
of study. 

Consider a diffusing particle placed initially at a position xq. The Fokker-Planck 
equation for the probability density P{-;t) of its position at time t is 



dt 



Pix;t) 



d_ 

dx 



d_ 

dx 



+ v{x) 



P{x;t) - &{x)P{x;t) 



where v and a are, respectively, the drift and the absorption (killing) rate of the 
diffusion. The transformation 

P{x;t) =e-ifo''^y^'^yilj{x;t) 

leads to the Schrodinger equation 

d2 



^ / 



HiP for := + ^ - - 
dx^ A 



When V and a are random, the expected value of P{xo; t) turns out to be indepen- 
dent of the starting point xq; it is related to the density of states N'{E) associated 
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with the disordered system via 

/>oo 

P{xo;t) = / dEN'{E)e'^^' . 

By standard Tauberian arguments, the long-time behaviour of this expected density 
is therefore completely determined by the low-energy behaviour of the density of 
states: 

lnA^(£') £'^as£'^0+ => lnP(xo;t) 1'' as t ^ oo . 

The Kotani case V — w corresponds to the choice v = and & = w. Here, 
the decay of the return probability density follows from the decay of the total 
probability due to absorption. For finite Levy measures, the Lifshits singularity 
(|6.23p leads to the well-known exponent = 1/3 [31]. On the other hand, when 
the Levy measure is infinite, the behaviour (|6.24p leads to the exponent 

1 + a 



The supersymmetric case V = w'^ /A — w' /2 corresponds to the choice v = w and 
a = 0. In this case, there is no absorption, the total probability is conserved, and 
the decay of the return probability density is due to the expected drift of the Levy 
process: 

/•oo 

E {W{x)) =x ym{dy) = 2ax . 
Jo 

For a pure deterministic drift, v ^ 1 and the return density decays exponentially. 
For a non-zero Levy measure, however, the fluctuations in the drift bring about a 
slower rate of decay. For example, in the case of a finite measure, we find, as in 
Kotani's case, v ~ 1/3; for the alpha-stable case, we find 

1 

3- 2a 

For the singular measure (|6.2ip . the spectral singularity (|6.22p leads to the prob- 
ability decay 

\nP{xo]t) ~ -t^'^hi''^t. 
In this marginal case, the singular nature of the measure is only refiected in the 
logarithmic correction. This behaviour holds for both the random drift case and 
the random absorption case. 

7. Conclusion 

In this paper, we have developed the supersymmetric version of the Frisch-Lloyd 
methodology for computing the density of states of a system with Levy disorder. 
The main novelty is that, in the supersymmetric case, the complex Lyapunov ex- 
ponent may be expressed in terms of the positive solution of a difference equation. 
When the Levy process is non-decreasing, the complex Lyapunov exponent has a 
continued fraction expansion whose coefficients have simple expressions in terms 
of the Levy exponent, and the problem of computing the density of states reduces 
to a Stieltjes moment problem. One pleasing outcome was the discovery of a new 
solvable case. Our efforts to adapt to the supersymmetric case Kotani's semiclas- 
sical analysis of the low-energy regime have also borne some fruit, although our 
approach only gives the leading term. This analysis has shown in particular that, 
for infinite Levy measures, the low-energy behaviour of the integrated density of 



28 



ALAIN COMTET, CHRISTOPHE TEXIER, AND YVES TOURIGNY 



states is controlled by the small jumps of the process. This result is in some sense 
the dual of that obtained in an earlier study by one of us, namely that it is the tail 
of the Levy measure that determines the high-energy regime [71 [5] . 
The paper raises a number of questions that deserve further study: 

(1) Much of the work presented here has been concerned with the subordinator 
case. From the point of view of localisation, however, the most interesting 
processes are those with zero means, such as Brownian motion. Hence we 
require a less restrictive theory. 

(2) As mentioned earlier, exponentials of Levy processes may be viewed as the 
zero-energy case of the supersymmetric model. There has been a great 
deal of recent work on such exponentials, including their close connection 
with self-similar processes |38j and a generalisation to complex exponentials 
|22| . Our own findings encourage us to seek interpretations of some of these 
results in terms of disordered systems. 

Appendix A. Some Levy processes 

We recall in this appendix some basic facts concerning the distributional prop- 
erties of two important classes of Levy processes. 

A.l. The alpha-stable process. We study the subordinator characterised by the 
infinite Levy measure 

a dy 
dy) = — rnr^' 0<a<l. 

The fact that the Levy exponent 

A{0) = -(-16')" 

has a simple power law indicates that the law is stable, namely, 

W{x) * = ' x'/'^Wil) 

or, in other terms, that the probability density function fw{' ]x) of the random 
variable W{x) satisfies 

/M/(u;x)=x-i/"/iv(a;-i/"u;l) 

where 

(A.l) /^(^;i)^ /^e'^-^'^)". 

This distribution is a particular case of a more general two-parameter family of 
Levy stable laws with probability density function 

C^p[u)= [ ^ e'^""l''l°(^+'''''sn(2)tan[7rQ/2]) _ 

Jr 27r 

The parameter /3 G [—1, +1] tunes the asymmetry of the law; /3 = corresponds to 
a symmetric law, such as the Cauchy Ci^o or the Gaussian £2,0 laws, while /3 = -f 1 
corresponds to subordinators, whose distributions have support in M+ [101116]. The 
precise relationship between fwi' ; 1) and Ca,/3 is 

fw{u ; 1) - [cos(7ra/2)]-i/" ([cos(7ra/2)]-i/"u) . 
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-tU — t°^ cos TTCt 



The distribution ()A.ip may be analysed as follows: for u < 0, the analyticity 
of the integrand in the lower complex plane implies f{u) = 0. For m > we 
can deform the contour of integration in order to follow the branch cut along the 
positive imaginary semi-axis; the result is 

f°° dt 

(A.2) fw{u;l)^ / — sin(i"sin7ra) e 

Jo I" 

and so we deduce an asymptotic power law tail similar to the Levy measure 

(A.3) fw{u]l)^-—^ — T""^"" asw^cxj. 

r(l - a) 

The behaviour of the density function (|A.1I) for small u may be studied by using 
a stationary phase approximation. Let us write 

/w(a1-i/";i) ^A^/^y ^e^^(") where (^(z) = iz - (iz)« . 

Since Re[iy9(z)] < for Imz < 0, we can deform the contour of integration in order 
to pass through the stationary point zq ~ — ia^/'-^""-' such that (fi'izo) = 0; this 
leads to 

1 2-c a ., 

(A.4) fwju;!)^ u exp [-AaU i-°J as u ^ 0+ 

where A„ = (1 - a)a"/'^^-"^ and B„ = (1 - a)a-^''^^-"\ 

For example, we can verify that these behaviours are the correct ones in the case 
a — 1/2, since £1/2. the density of the so-called Levy distribution, admits the 
simple expression 

fw{u ; 1) = ^ u-3/2e-i/(4«) 1^^^^^ (y) . 

A.2. The gamma process. The gamma process is characterised by the infinite 
Levy measure [1] 

m(dy) = '^e-^y . 

Using 

we obtain the Levy exponent 

(A.5) A(6i) = -ln(l-i6'/5) . 

From the point of view of the statistical properties of the frequent small jumps, the 
gamma process can be understood as the a — > limit of the alpha-stable process. 
The distribution can be studied by considering the Fourier transform 

dz e^' 

fw{u]x) = 
A contour deformation leads to 

(A.6) /^(^,;2;) = -^«-ie"''"l(o,oo)H 

In particular, the moments are finite, given by 

T{n + x) 



27r (1 -|-iz/6)=' 



nwixY^ r(.) 
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